% -*- mode: octave -*- (octave mode for emacs)
function y=normeH1(t,u,du)
  z=u.^2+(gradient(u,1/(length(t)-1))).^2;
  %% z=2*(u.^2);
  y=sqrt(trapz(t,z));
  %%end
